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It is shown that, in the self-consistent quantum statistical Hartree-Fock approxi- 
mation, the number of electronic states localized on one nucleus is finite. This result 
is obtained on the basis of the general electron-nuclear model of matter and pro- 
vides convergence of the atomic statistical sum and finiteness of the "atom" size. 
In general approach the characteristic size of the "atom" is a function of density 
and temperature. However, it is shown, that in a wide range of thermodynamic pa- 
rameters, for relatively low temperatures, characteristic orbits and electron energy 
eigenvalues are independent of density and temperature. In this case, the sizes of 
the orbits are of order of the Bohr radius which is a minimal characteristic size in 
the system for typical parameters of plasma with atomic states. 
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I. INTRODUCTION 

Real matter (gases, liquids, solids, plasma) is constructed from nuclei and electrons in- 
teracting with each other according to the Coulomb law [1]. It would be very desirable, 
based on the known nucleus charge as an initial problem parameter and, based only on 
the Coulomb law, to construct a scheme for calculating macroscopic properties of different 
media. The consistent approach to implement such a scheme is based on the methods of the 
diagram technique of the perturbation theory [1,2]. The main difficulties in the realization of 
this approach are caused by the necessity of the simultaneous consideration of both " local- 
ized" states of electrons and nuclei ("atoms", "molecules", etc.), whose description requires 
accurate summation of rows of the perturbation theory, and " delocalized" states which are 
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quite successfully described within the lowest orders of the perturbation theory [1,2]. As a 
result, "localized" states for weakly nonideal plasma are considered within the second virial 
coefficient [1-4]. In this case, there arises the known problem of the convergence of the sta- 
tistical sum of the hydrogen atom, which is directly related to its size [2]. A large number 
of works are devoted to the solution of this problem, where various procedures and physi- 
cal mechanisms of "cutoff' of the statistical sum of the atom were proposed, including the 
papers [5-12] which have already become classical (see [1-3] and references therein for more 
details). In this case, the following circumstance is very surprising. As is known [13,14], 
the self-consistent Hartree-Fock approximation and its various generalizations are basic in 
calculations of the structure and energy spectra of atoms and molecules in quantum chem- 
istry. Despite the known statement that the self-consistent Hartree-Fock approximation is 
the best single-particle approximation in quantum statistics [15], it is almost not used in the 
study of thermodynamic properties of Coulomb systems, where "localized" states should be 
taken into account [1]. This paper is devoted to the solution of the problem of convergence of 
the statistical sum of the hydrogen atom and its size finiteness, based on the self-consistent 
quantum-statistical Hartree-Fock approximation for the Coulomb system. 



II. TWO-COMPONENT ELECTRON-NUCLEAR PLASMA: ADIABATIC 
APPROXIMATION AND CLASSIFICATION OF THE ELECTRON STATES 

We consider the two-component Coulomb system (CS) in volume V at temperature 
T within the canonical ensemble consisting of electrons (subscript e) and nuclei (sub- 
script c), taking into account the quasineutrality condition ^a=e c^a.e''^a'^ = 0. Here 
ni°^ = {NaY^'^yV is the average density of the number of particles; Na is the operator 
of the total number of particles of type a, characterized by mass nia, the average density 
of the number of particles Ua, and charge ZaC; angle brackets (. . . mean averaging over 
the canonical ensemble. For such a CS, the Hamiltonian H^^ — H^e + i?cc + -f^ec in the 
secondary quantization representation is written as [1] 

^faa = / V'a (r) V'Va(r)rfr + Ka, H^c = j u^d^i - Fa) A^e(ri) A^c(r2)c?ridr2, (1) 

Vaa^ ^ / Uaa{ri - r2)'0a (rijV'a (r2)'0a(r2)'0a(ri)(iridr2, Uaa{r) = °^ , (2) 
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Here Uaai^) is the potential of the Coulomb interaction of particles of types a and b, 
■0a (r) and ■0a(r) are the field creation and annihilation operators, Na{r) — ■0+(r)'0a(r) is the 
density operator of the number of particles of type a, which, in the coordinate representation, 
is written as Na{r) = '^f=iS{r~'R,f). The vahics arc the coordinates of particles of 
type a. We emphasize that the total number of particles Na of each type in the canonical 
ensemble is specified and is the C-number, Na — J Na{r)dr — (Na)^'^'^^. The free energy 
F*^^ of the considered quasi-neutral CS consisting of the set number of electrons and 
nuclei Nc in volume V at temperature T, F^^ — —TlnSp |exp ^— -^f-j |, in the adiabatic 
approximation for nuclei, is written as 

FOS ^ _TlnSp^ jexp } , Fe. = -Tln^p, jexp (^-^hl±^^ | . (3) 

The quantity F^.^ (3) is the free energy of the subsystem of electrons being in the external 
field of stationary nuclei. In this case, the operator of the density of the number of nuclei 
nc{v) in the Hamiltonian H^c (1) is the C-number. At this stage, when considering the sub- 
system of electrons in the external field of stationary nuclei, we proceed from the equivalence 
of the large canonical and canonical ensembles, i.e., we suppose that the equality [16] 

Fee = ^ec + Ate ( A^e) , ^ec = -T In Spe | CXp 

for the free energy Fee (3) is valid. In (4) the function flee is the thermodynamic potential of 
the subsystem of electrons with the chemical potential /le in the external field of stationary 
nuclei, whose value depends on the nucleus coordinates {R^} as the external parameters, 
(Ai'e)^*^^^ is the average number of electrons as a function of temperature T, volume V, and 
chemical potential fie- In this case, the quasineutrality condition takes the form (Ne)^^^'' — 
ZcNc and is used to determine the chemical potential fie of electrons. 

Then, to describe the field operators "0^ (r) and ipa{^) in the representation of occupation 
numbers [13], it is required to choose a complete set of wave functions (r|{(5(7}, {R?}) 
characterizing single-particle electronic states with the energy E{{Q}) 

{Qa} {Qa} 

in the external field of nuclei. Here a'^q^y and a^Qa} are the creation and annihilation 
operators of electrons in the state (r|{(5cr}, {R?}), {Qa} is the set of quantum numbers, 
including the spin number a (hereafter, for simplicity, spin indices are not considered). In 
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the problem under consideration the set of coordinates {R^} of nuclei plays the role of 
eigenvalue parameters. In this case, single-particle electronic states in the external field of 
nuclei are classified into two groups [17]. One (A) is localized (subscript "loc") states char- 
acterized by exponential spatial decrease with distance from localization centers. Among 
them are the single-center bound states of electrons in atoms, two-center ones in molecules, 
etc. Another group (B) is delocalized (subscript "del") states in which the electron prop- 
agates throughout the system. The typical example of delocalized states are plane waves 
(r|q), (r|q) = Ij^fV exp(iqr), e(g) = |^ which describe the behavior of the free electron 
with momentum and energy e(g). Thus, in the optimum case, the complete set of wave 
functions for the representation (5) should describe both localized and delocalized states of 
electrons. However, there are many problems in this way. The first problem is "physical"; it 
is associated with the fact that localized electronic states vary as thermodynamic parameters 
of the CS are varied; this leads, in particular, to the Mott effect [18], i.e., disappearance 
of localized states. This means that the quantum-mechanical consideration is insufficient 
for determination of localized states; quantum statistics effects should also be considered. 
The second problem, to the solution of which this paper is devoted, is associated with the 
necessarily finite number of electronic states localized on one center. By the example of the 
hydrogen atom considered within quantum mechanics ignoring quantum statistics effects 
[13], it is known that the number of bound electronic state is infinite, which leads, in turn, 
to the divergence of not only the statistical sum for an individual atom, but also to the 
unlimited atomic size at any nonzero temperature [2]. Based on combinatory reasons, E. 
Fermi [7] has shown that the introduction of the finite size of atoms results in an exponential 
limitation of the statistical sum. It is easy to see that the statistical sum limitation to any 
maximum value of the principal quantum number n^nax is also a simultaneous introduction 
of the finite size of the atom. The absence of limitation, i.e, the case n — oo, is the case 
corresponding to an infinite-size atom, which is meaningless (the size of the n-th orbit of 
the hydrogen atom is it!„ = n^oo, where is the Bohr radius) [19]. Hence, it is impossible 
to correctly use the "single-center" (in the sense of the localization center) approximation 
in the description of localized states. In this case, the transition to the "single-center" 
approximation is meaningful when the condition 



< Ra >« r > 
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is satisfied, where < Ra > is the "atom" size and < r >c is the average distance between 
nuclei (see below). 

The consideration of the quantum statistics effects (with the opportunity to pass to the 
thermodynamic limit) means that localized electronic states for a separated nucleus should 
be considered taking into account that the CS, as a system of many identical particles, 
contains many other electrons (electronic states). However, the multiplicity of noninter- 
acting electrons is already formally taken into account in the used formalism of secondary 
quantization (1), (2), (5). 

The divergence directly follows from the definition of the electron distribution function 
fe{E{{Q})) in states {r\{Q}, {K^}) with energy E{{Q}) 

ME) = {al^^HQ^r^^ = jexp + ij" , (iV.)(«^) = J^UE), (7) 

^ ^ ^ J {Q} 

when using energy levels of the hydrogen atom as energies. It remains to hope that the "ef- 
fective" consideration of the electron-electron interaction while retaining the single-particle 
consideration will make it possible to solve the problem of the finite number of bound states 
in the "atom", i.e., in the "single-center" approximation for wave functions (it is clear that 
the "multicenter" problem for localized electronic states is in principle unsolvable). 



{v\{loc}, {RO) ^ (r - R^|n) = (^„(r - R^), E{n, {loc}) ^ (8) 

In this case, the "atom" size in the "single-center" approximation, taking into account 
(7) and (8), can be defined for quantum statistics as follows 

< Ra f r'J2fe{En)\<Pn{r)\'dr. (9) 

We emphasize that such quasiparticles - "atoms" are not identical to atoms commonly 
used in quantum mechanics. The difference is caused by the fact that any of one of a 
very large (infinite in the thermodynamic limit) number of electrons can be found in the 
corresponding electronic state of the "atom" due to the identity, which is reflected in the 
statistical description of such a system. Thus, electronic states in " atoms" will depend on 
the thermodynamic parameters, i.e., the temperature T and average density of the number 
of electrons nf^ (or on the chemical potential //g) [2]. 
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We note that electronic states are considered in the random field of nuclei. In this case, 
the role of quantum numbers for both "localized" (with discrete energy spectrum) and 
" delocalized" (with continuous energy spectrum) electronic states is played by the energy 
itself of the corresponding state and the spin quantum number [20]. For "localized" states, 
this set of quantum numbers can be complemented by coordinates of localization centers 
(nuclei) {R^} [20], since electronic states with this energy can become localized at different 
points of the system. In this sense, we can speak of degeneracy with respect to the quantities 
{R^} which, therefore, can be considered as "quantum" numbers. In this case, the term 
"quantum number" should be understood in a certain conditional meaning: the vectors {R^} 
and their components are not eigenvalues of any operator having an effect on "electronic" 
variables. 

III. FINITE ATOM IN HARTREE-FOCK APPROXIMATION 

Further, let us consider the system of electrons in the external field of nuclei within the 
self-consistent Hartree-Fock approximation for quantum statistics [15]. Then the complete 
set of wave functions (r|{(5}, {Rf}) and the corresponding energy spectrum E{{Q}) of the 
multicenter problem are determined from the set of equations 



In (10), the term with function (ne(r)) = (ne(r, ri)) |r=ri corresponds to the Hartree 
approximation. The term with function (ne(r, ri)) = Yl{Q}{^\{Q}){{Q}\^i)fe{E{{Q})) in 
(10) corresponds to the Fock approximation. Nonlocal structure of the Fock term in (10) 
significantly complicates the problem. However, its consideration is conceptually necessary. 
First, in the Hartree approximation, the states (r|{Q},{R^}) are not orthogonal to each 
other [13]. Second, the consideration of the Fock approximation allows elimination of the 
so-called "self-action". As it is easy to check from (10), the state (r|{(5}) in the Hartree-Fock 
approximation is determined only by other states, in contrast to the Hartree approximation 
[13]. The set of equations (10) includes not only localized, but also delocahzed states, to 
describe which the "single-center" approximation is inapplicable. Taking into account that 
we are interested in justification of the finiteness of the number of localized states, let us 
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simplify the set of equations (10) under the assumption that the number of delocahzed 
states is very small. This can be done using the properties of the Fermi-Dirac distribution 
(7) for the function fdE) in (10). Let us perform the further consideration without loss of 
generality immediately for the case of hydrogen, Zc = 1, {N^)^'^^^ = Nc- Let us suppose that 
for the energy Eq of the ground single-particle localized state, energies Egxc of each of other 
states, including delocalized states, and the chemical potential /ig of electrons, the following 
relations are valid 

£;0<0, |^o|>T, |^o-^exc|>r, Eo<l^e<0- (H) 

Then, according to (7), 

MEo) - 1 - exp - 1, /e(£;e.c) = exp ^^^^^ « 1. (12) 

We further take into account that the number of single-particle locahzed states is pro- 
portional to the number of localization centers; in the "single-center" approximation, it is 
proportional to the number of nuclei Nc- We also assume that the delocalized states are 
close to corresponding plane waves. Then, with required accuracy, we can consider that 
the energy levels of delocalized states are defined by the energy e{q) of the free electron, 
E{{del}) e{q). Thus, taking into account 

(iV,)(«^) ^Y.feiE{{Q})) ^ N,J2fe{En) + $]/e(6(g)), (13) 

{Q} " 1 

and (7)-(12), wc find the equation for determining the chemical potential /j.^ of electrons in 
the "single-center" approximation for localized states, 

w)'-' = {i - exp (^) + x:-p (^) } + E«p (s^) • 

Then, taking into account (11) and the equality 

(S^) - V/ {^) A ^ (S) 

from (14), we find 
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Thus, conditions (11) for the chemical potential are satisfied when 

E-pf^)+^exp(f)<l. (17) 



Taking into account (11), inequahty (17) is valid for a hmited number of localized states. 
According to (11), (12), (14), the contribution of delocalized states in the self-consistent 
Hartree-Fock approximation (10) can be neglected in a very wide range of thermodynamic 
parameters satisfying the condition X^wt^o^-^P '^TJ^' When determining localized 

states and the corresponding energy spectrum, we can consider with required accuracy that 
the ground state is "populated", fe{Eo) = 1 and all excited states are "free", fe{Eexc) — 0. 
Thus, when relations (11), (16), (17) are valid, the set of equations of the self-consistent 
Hartree-Fock approximation (10) can be written for localized states in the "single-center" 
approximation (8) (the nucleus with charge Zc — 1 is at the coordinate origin), 

2me 

Uee{r - ri)cpo{r)cpl{ri)cpn{ri)dri = E„cpn{r), (18) 



+ J itee(r - ri)|(/7o(r)|^rfri| (/7„(r) - 



It immediately follows from (18) that the wave function (po{^) and the ground state energy 
are identical to the similar quantities for the electron ground state in the hydrogen atom 
[13], (po{r) = (ttOo)""^ exp (^-^) , Eq = where oo is the Bohr radius. The "excited" 

locahzed states {n > 1) are characterized by the energies En and differ from corresponding 
states in the hydrogen atom, which makes it possible, in particular, to remove the known 
"Coulomb degeneracy" of energy levels [13]. Taking into account the explicit form of the 
Fourier components of the Coulomb interaction potential (2), we can calculate the Hartree 
term, i.e., the last braced term in (18), 

UH{r) = f Uee{r - ri)|(^o(r)|'(iri ^ - - - (l + —] ei^p ( -—] . (19) 

Further analytical calculations seem impossible. As in quantum-mechanical calculations, 
in view of the nonlocality of the Fock term, equations (18), taking into account (19), are 
solved using iterations [13,14]. However, by analogy with the states in the hydrogen atom 
[13], we can assume that the limit relations 

(a T \ 
— , < Qn <1. (20) 
«0 / 
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are valid for functions (fin{r){n > 1). 

In this case, the Fock term in (18) and (19) can be written in the "local" form 

J Uee{r - ri)(po{r)ipl{ri)(pn{ri)dri = u'^{r)(pn{r), (21) 
Upir) = Mr)ip;^\r) J Mee(g) exp(iqr)0on(g)^^^, (22) 

0on(?) = J Vo{r)Vn{r) exp(-iqr)(ir. (23) 

According to (20) and (23), 

lim (j)Qn{q) = 0, lim0on(g) = (pon{q = 0) = / (po{r)(pn{r)dr = 0. (24) 

due to the orthogonality of the functions (^o(^) and 93„(r). Then it immediately follows from 
(21)-(24) that 

— — ] , 7n > 0. (25) 
oo / 

Hence, we come to the conclusion that the function (/?n(r) is defined as the wave function 
satisfying the Schrodinger equation. 



2m. 



e 



C/„(r^O) = --, C/„(r^oo)~expf-^y (26) 
r V ^0 / 

Thus, the initial assumptions (20) are valid. Furthermore, the asymptotic behavior of the 
potential Un{r) is similar to the behavior of the known Yukawa potential [21]. As is known, 
at certain parameters of the Yukawa potential, it will not contain the bound (localized) 
state [22]. As applied to the problem under consideration, this means that the number of 
locahzed states in the hydrogen "atom" is limited. 

Using the results obtained above, it is easy to see that the free energy F'"^ (3) of the 
Coulomb system under consideration equals 



nC5 



-NcT < 1 + In 



, /e(T) = exp ( ) (27) 



10 



if the conditions (6), (11), (16), (17) are valid. 

The relation (27) corresponds the free energy of the ideal gas of new quasiparticles - 
"atoms", in which electron is in the ground state (see, e.g., [16]). 

Summarizing the above consideration, it can be assert that the statistical sum and the 
hydrogen "atom" size are finite when using the self-consistent quantum-statistical Hartree- 
Fock approximation. In this case, in the above approximation (11) for the chemical potential 
of electrons, the "atom" size (9) is independent of the density of the number of particles 
and of temperature in the system and is defined only by the Bohr radius in a wide range 
of thermodynamic parameters. For low temperature the developed approach provides the 
rigorous derivation of the atomic free energy (27), written earlier on the basis of physical 
reasons. On the basis of the above derivation, it is evident that the finiteness of the number 
of localized states leads to the finite free energy of the CS at arbitrary temperature. 

Acknowledgments 

The authors thank W. Ebeling and A.L. Khomkin for the useful discussions. This study 
was supported by the Netherlands Organization for Scientific Research (NWO), project no. 
047.017.2006.007, and the Russian Foundation for Basic Research, projects no. 07-02-01464- 
a and no. 10-02-90418-Ukr-a. 



[1] W.-D. Kraeft, D. Kremp, W. Ebeling and G. Ropke, Quantum Statistics of Charged Particle 

Systems (Plenum Press, New- York, 1986) 

[2] W. Ebeling, W.-D. Kraeft and D. Kremp, Theory of Bound States and Ionization Equilib- 
rium in Plasmas and Solids (Akademie-Verlag, Berlin, 1976) 
[3] W. Ebeling, Ann. Physik, Leipzig 22, 383, 392 (1969) 
[4] A.N. Starostin and V.C. Roerich, JETP 100, 165 (2005) 
[5] K.F. Hcrzfcld, Ann.der Phys. 19, 380 (1916) 
[6] M. Planck, Ann.der Phys. 75, 673 (1924) 

[7] E. Fermi, Z.Phys. 26, 54 (1924), Nuovo Cimento 11, 157 (1934) 

[8] L. Brillouin, Les Statistiques Quantiques et leurs Applications (Les Press.Univers.Prance, 



11 



Paris, 1930) 

[9] K.H. Riewe and R. Rompe, Ann.der Phys. Ill, 79 (1938) 
[10] A. Unsold, Z.Astrophys. 24, 79 (1948) 
[11] D.R. Iglis and E.Teller, Astrophys. J. 24, 355 (1948) 
[12] A.I. Larkin, JETP 38, 1896 (1960) 

[13] A. Messiah, Quantum Mechanics (North-Holland, Amsterdam, 1961) 

[14] C. Froese Fisher, The Hartree-Fock Methods for Atoms (Wiley-Interscience, New- York, 

1977) 

[15] D.Kirzhnits, Field Theoretical Methods in Many-body Systems (Pergamon, New York, 
1967) 

[16] L.D. Landau and E.M. Lifshitz, Statistical Physics, part 1 (Butterworth-Heinemann, Ox- 
ford, 1980) 

[17] V.B. Bobrov and S.A. Trigger, Solid State Comm. 56, 29 (1985) 

[18] N.F. Mott, Metal-Insulator Transitions (Taylor and Francis, London, 1974) 

[19] L.P. Kudrin, Statistical plasma physics (Atomizdat, Moscow, 1974) 

[20] V.L.Bonch-Bruevich, R.Enderlein, B.Esser, R.Keiper, A.G.Mironov and I.P.Zvyagin, Elec- 
tronentheorie ungeordneter Halbleiter (VEB Deutch.Verlag, Berlin, 1984) 

[21] H. Yukawa, Proc. Phys. Math. Sac. Japan 17, 48 (1935) 

[22] S. Flugge, Practical Quantum Mechanics (Springer- Verlag, Berlin, 1971) 



